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ABSTRACT 

We present the results of numerical simulations which show the formation of blue 
compact dwarf (BCD) galaxies from merging between very gas-rich dwarfs with ex- 
tended HI gas disks. We show that dwarf-dwarf merging can trigger central starbursts 
and form massive compact cores dominated by young stellar populations. We also 
show that the pre-existing old stellar components in merger precursor dwarfs can be- 
come diffuse low surface brightness components after merging. The compact cores 
dominated by younger stellar populations and embedded in more diffusely distributed 
older ones can be morphologically classified as BCDs. Since new stars can be formed 
from gas transferred from the outer part of the extended gas disks of merger precursors, 
new stars can be very metal-poor ([Fe/H] < — 1). Owing to very high gaseous pressure 
exceeding 10^ fc^ (where ks is the Boltzmann constant) during merging, compact star 
clusters can be formed in forming BCDs. The BCDs formed from merging can still 
have extended HI gas disks surrounding their blue compact cores. We discuss whether 
tidal interaction of gas-rich dwarfs without merging can also form BCDs. 
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1 INTRODUCTION 

. Blue compact dwarf galaxies (BCDs) are defined as dwarfs 
■ having Mb ^ —18 mag, optical sizes smaller than 1 kpc, 
and spectra similar to HII regions of spiral galaxies (Thuan 
& Martin 1981) and their formation and evolution processes 
have long been discussed by many authors. Some gas-rich 
BCDs dominated by very metal-poor, young stellar popula- 
tions (e.g., I Zw 18) were once suggested to be undergoing 
their first starbursts (e.g., Searle et al. 1973; Aloisi et al. 
1999; Thuan et al. 1999), recent observations have however 
confirmed that most of BCDs have underlying older stel- 
lar populations at least a few 10^ yr old (e.g., James 1994; 
Cairos et al. 2001; Amon'n et al. 2007). Although optical 
structures (Gil de Paz et al. 2003), HI properties (van Zee et 
al. 1998), and chemical abundances (e.g.. Hunter & Hoffman 
1999) in BCDs have been discussed in terms of evolutionary 
links between BCDs and other types of dwarfs such as dwarf 
spheroidal (dSphs), ellipticals (dEs), and irregulars (dirrs), 
physical connections between these dwarfs remain unclear. 

A growing number of observations have recently sug- 
gested that the origin of BCDs can be closely associated 
with galaxy merging between low- mass dwarfs (e.g., Noeske 
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et al. 2001; Ostlin et al. 2001). For example, Pustilnik et al. 
(2001) investigated environments of BCDs (e.g., the pres- 
ence or the absence of nearby companion galaxies) and sug- 
gested that at least 80% of BCDs in the studied sample 
have evidence of their starburst components being triggered 
either by tidal interaction or by galaxy merging. Further- 
more, observational studies of the six most metal-deficient 
starbursting dwarfs (most of which are BCDs) by Giant Me- 
trewave Radio Telescope (GMRT) have reported that they 
show clear signs of interaction and merging in their HI prop- 
erties (Chengalur et al. 2008). However, no theoretical and 
numerical studies have yet provided the detailed predictions 
on physical properties of BCDs formed from "dwarf-dwarf 
merging" . Therefore, it remain unclear whether the observed 
properties of the BCDs can be consistent with the merger 
scenario of BCD formation. 



The purpose of this paper is thus to present the results 
of numerical simulations which show physical properties of 
BCDs formed from dwarf-dwarf merging. Based on the nu- 
merical results, we discuss (i) the origin of young, metal-poor 
stellar populations in BCDs, (ii) possible evolutionary links 
between BCDs, dSphs, dEs, and dIrrs, and (iii) the roles of 
tidal interaction in the formation of BCDs. 
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Table 1. Model parameters for N-body simulations 
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disk 
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M2 


disk 


1.0 


5.00 


0.1 


1.0 


M3 


spheroid 


1.0 


5.00 


0.1 


1.0 


M4 


disk 


0.4 


1.25 


0.04 


4.0 



Orbit types f 


m2 


„ h 





Comments 


PR ~ 


0.5 


1.0 


1.0 


merging 


PP 


1.0 


1.1 


10.0 


interaction 


PP 


1.0 


1.1 


10.0 


interaction 


PR 


0.5 


1.0 


1.0 


gas-poor 



" The initial morphology of old stars: "disk" and "spheroid" represent the dwarf disk and spheroidal models, respectively. 
^ The initial total stellar mass of a dwarf in units of 10^ Mq . 

The initial stellar size of a dwarf in units of kpc. 

The initial gas mass in units of 10^ Mq. 
^ The initial size ratio of gas disk to stellar one. 

f "PR" and "PP' describe prograde-retrograde and prograde-prograde orbits, respectively. 

9 The mass ratio of interacting or merging dwarfs. 

^ The orbital eccentricity of interacting or merging dwarfs. 

* The pericenter distance of interacting or merging dwarfs in units of kpc. 



2 MODEL 

We investigate chemodynamical evolution of gas-rich merg- 
ers between dwarfs embedded in massive dark matter ha- 
los by using our TREESPH chemodynamical codes that 
combine key ingredients used in our previous codes: chemi- 
cal evolution and supernovae feedback effects (Bekki & Sh- 
ioya 1999), hydrodynamics for strongly self-gravitating gas 
(Bekki 1997), and formation of globular clusters (GCs) and 
field stars (Bekki et al. 2002) . This code enables us to investi- 
gate chemical and dynamical evolution of forming BCDs in a 
self-consistent manner. We adopt the Burkert profile (Burk- 
ert 1995) for the radial density profile of the dark matter 
halo of a dwarf, because it can be consistent with rotation 
curve profiles of dwarfs (Burkert 1995). The total mass of the 
dark matter (Mdm) for a dwarf is set to be 8.0 x lO^M© for 
all models and the dark matter halo has a large core radius 
(rc) of 3.5 kpc and the truncation radius of 3.4rc (Burkert 
1995). The old stars in the dwarf are assumed to have ei- 
ther a disky distribution (referred to as "disk model") or a 
spheroidal one ("spheroidal model"). 

The stellar component of a dwarf in the disk model is 
described as a purely disk with the initial mass of Ms and 
the disk size of Ts- The radial (R) and vertical (Z) density 
profile of the initially thin disk are assumed to be propor- 
tional to exp{—R/Ro) with scale length Ro — 1 kpc and to 
sech^(Z/^o) with scale length Zq = 0.2Ro, respectively. For 
the spheroidal models, the stellar spheroids are assumed to 
have the projected radial density profiles exactly the same as 
those for the stellar components in the above disk models. 
In both disk and spheroidal models, a dwarf has a thin gas 
disk and the gaseous disk is assumed to have the same scale 
length as that of the stellar one. The size- and mass-ratios of 
the gaseous disk to the stellar one are set to be free param- 
eters and described as Sg and /g, respectively: the size (rg) 
and mass (Mg) of the gas disk are SgVs and fgMs, respec- 
tively. The HI diameters of gas-rich galaxies are generally 
observed to be larger than their optical disks (Broeils & van 
Woerden 1994). A small fraction of low luminosity galaxies 
have HI gas envelopes extending out to 4-7 Vs (e.g.. Hunter 
1997) with HI mass to light ratios up to ~ 2OM0 Lg^ (War- 
ren et al. 2004). Guided by these observations, we investigate 
models with 1 ^ Sg 4 and 0.1 ^ /g ^5 2. 

Star formation is modeled by converting the collisional 



gas particles into coUisionless new stellar particles accord- 
ing to the algorithm of star formation described below. We 
adopt the Schmidt law (Schmidt 1959) with exponent 7 — 
1.5 (1.0 < 7 < 2.0, Kennicutt 1998) as the controlling pa- 
rameter of the rate of star formation. The stars formed from 
gas are called "new stars" whereas stars initially within a 
disk are called "old stars" throughout this paper. Chemical 
enrichment through star formation and supernova feedback 
is assumed to proceed both locally and instantaneously in 
the present study. The values of chemical yield and return 
parameter are 0.002 and 0.3, respectively, and the initial 
gaseous metallicity is [Fe/H] — —1.6. 

Nearby BCDs (e.g., NGC 1705) are observed to have 
compact young star clusters (SCs) that can finally evolve 
into GCs (e.g., Meurer 1993). We try to investigate whether 
BCDs formed from dwarf-dwarf merging and interaction 
can have SCs. We adopt a plausible assumption that if gas 
pressure (Pgas) can exceed a threshold pressure (Pth), such 
high pressure of ISM can induce the global collapse of giant 
molecular clouds to form massive compact star clusters cor- 
responding to SCs (Jog & Solomon 1992; Elmegreen & Efre- 
mov 1997) . We show the results of the models with Pth = 10^ 
fca, where fce is Boltzmann's constant. The present results 
do not depend so strongly on Pth for Pth = lO'^-lO'^fca. New 
stars formed from gas with Pgas > 10"" fee are identified as 
SCs rather than field stars in the present study. We mainly 
investigate mass fractions of SCs among all new stars for 
each model. 

The mass ratio of the two merging or interacting dwarfs 
(7712), the pericenter distance (^p), and the eccentricity (ep) 
are assumed to be free parameters. The orbit of the two 
dwarfs is set to be the xy plane and the distance between 
the center of mass of the two dirrs is 20 kpc. The spin of 
each galaxy in a merger is specified by two angles Oi and 
where suffix i is used to identify each galaxy. Oi is the angle 
between the z axis and the vector of the angular momentum 
of a disk. 0^ is the azimuthal angle measured from the x axis 
to the projection of the angular momentum vector of a disk 
onto the xy plane. 

Although we run many models with different param- 
eters with different m2, /g, Sg, and orbits of mergers, we 
show only four representative models (Ml, M2, M3, and 
M4) with two different orbital configurations: the highly in- 
clined prograde-retrograde ("PR") orbit with rp — 1.0 kpc. 
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Figure 1. Distributions of old stars (magenta), gas (cyan), and 
new stars (yellow) projected onto the x-z plane at T = 0.50 Gyr 
(left) and T = 1.33 Gyr (right) in the model Ml. The time T 
shown in the upper left corner for each panel represents the time 
that has elapsed since the simulation starts. The numbers shown 
in the upper right corner for each panel represents the star for- 
mation rate (Mq yr~^). The bar in the lower left corner for each 
panel measures 2 kpc. 
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Figure 2. The final projected radial density profiles of the simu- 
lated BCD for old stars (thin) and new ones (thick) in the model 
Ml. The profiles of the merger remnant well dynamically relaxed 
are shown here. 



ep = 1.0, 01 = 30°, 02 = 120°, 4>i = 90°, 4>2 = 30° and the 
prograde-prograde ("PP") one with rp = 10.0 kpc, ep = 1.1, 
01 = 0°, 6I2 = 0°, 01 = 0°, (^2 = 0°. We show the results of 
the model Ml in detail, because it shows typical behaviors 
of BCD formation from dwarf-dwarf merging in the present 
simulation. We also show the results of the model M2 and 
M3 in which two dwarfs are tidally interacting with each 
other without merging, because this comparison of the two 
models enables us to clarify how interacting galaxies can be- 
come BCDs. The parameter values for the four models are 
shown in the Table 1. The total particle number used in a 
major merger model is 220000 for collisionless particles and 
30000 for collisional ones. The fixed gravitational softening 
length for dark matter, old and new stars, and gas are set 
to be 0.35 kpc, 0.04 kpc, and 0.135 kpc, respectively. We 
confirm that stellar and gaseous disks are stable in isolated 
models (i.e., with no merging and tidal interaction) even if 
gas mass fraction are quite large: this is due to the very 
extended gas disks in the present models. 
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Figure 3. The final radial profile of mean metallicities for new 
stars within the central 5kpc of the merger remnant in the model 
Ml. The la dispersion in the metallicities for each radial bin is 
shown by an error bar. The mean metallicity (^m) of new stars 
is [Fe/H] = —0.93 for this model and Zm depends totally on the 
adopted initial gaseous metallicities and chemical yield parame- 
ters in the present models. 

3 RESULTS 

Fig. 1 shows mass distributions of different components (old 
stars, gas, and new stars) in the merging model Ml at two 
different epochs with star formation rate (SFR) of 1.O5M0 
yr~^ and O.OSM© yr~^. SFR reaches its maximum value of 
4.I3M0 yr~^ and about 80% of initial gas mass is consumed 
within 1.3 Gyr. In the late stage of merging between two gas- 
rich dwarfs with extended gas disks (T = 0.50 Gyr), nuclear 
starbursts can be triggered owing to efficient inward trans- 
fer of gas from the outer HI disks of the merger progenitor 
dwarfs. As a result of nuclear starbursts, the compact core 
dominated by young stars can be formed and embedded in 
old stars and gas. 

The spatial distribution of old stars clearly shows a dis- 
turbed morphology in the outer part of the merger {T = 0.50 
Gyr). If the stellar population synthesis models for metal- 
licities of [Fe/H] = —1.28 and ages of 0.5 Gyr by Vazdekis 
et al. (1996) are applied to the central core dominated by 
young stars, then B — R is estimated to be ~ 0.6, which is 
consistent with the observed range of B — in BCDs (Gil De 
Paz et al. 2003). This merger with a blue compact core dom- 
inated by young stars and an outer disturbed morphology 
can be thus identified as a BCD. 

Fig. 1 clearly shows that after strong starbursts, the 
dwarf-dwarf merger can soon get dynamically relaxed to 
form a flattened spheroid with a central compact core and a 
HI envelope (T — 1.33 Gyr). The HI envelope with the mass 
of 1.6 X lO^M© (i.e., M^/M^ = 0.27) has a disky distribution 
with its outer part strongly warped. The SFR after strong 
starbursts can still become relatively high (e.g., ~ O.IM© 
yr~^ at r = 1.2 Gyr) in a sporadic way owing to the pres- 
ence of the gas disk. This merger remnant with a compact 
core, a more regular morphology, and a higher SFR therefore 
can be also identified as a BCD. 

The projected radial density profiles of old and new 
stars shown in Fig. 2 clearly demonstrate that the merger 
remnant is dominated by new stars formed from nuclear 
starbursts in the central 1 kpc. If we adopt Mb /Lb (i.e., ratio 
of stellar mass to light in B-band) of 1.0, B-band surface 
brightness (/ib) distribution of old (new) stars ranges from 
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Figure 4. The time evolution of star formation rate (SFR) for 
the tidal interaction models M2 (solid) and M3 (dotted). 



- 24 (~ 22) mag arcsec"^ at - 0.1 kpc to ~ 29 (~ 29) 
mag arcsec"^ at i? ~ 2 kpc in this model. Since Me /Lb 
is significantly smaller than 1 for stars younger than 1 Gyr 
(e.g., Vazdekis et al. 1996), hb in the central region of the 
simulated BCD would be significantly higher than the above. 

Fig. 3 shows that metallicities in [Fe/H] for new stars 
within the central 5 kpc of the merger remnant ranges from 
~ —1.3 to ~ —0.8: the simulated BCD has a significant neg- 
ative metallicity gradient of new stellar populations. Owing 
to inhomogeneous chemical mixing during starbursts, the 
metallicity dispersion at a give radius becomes significantly 
large (~ 0.2 dex). The derived very small metallicities are 
due essentially to the inward transfer of metal-poor gas ini- 
tially in the extended HI disks of merger progenitor dwarfs. 
It should be however stressed here that the metallicities of 
new stars depend strongly on the adopted initial gaseous 
metallicities in the present models. These results imply that 
younger stellar populations in BCDs, which consists of their 
blue compact cores, can be very metal-poor ([Fe/H] < —1), 
because gas that formed the younger populations originate 
from chemically pristine HI gas initially in the outer part of 
merger progenitor dwarfs. 

Fig. 4 shows time evolution of SFR for the tidal interac- 
tion models M2 and M3 in which the initial dwarf is modeled 
as a disk and a spheroid, respectively. As shown in Fig. 4, 
a secondary starburst can be triggered to form a compact 
core dominated by young stars only in the disk model (M2): 
nuclear starbursts that form blue compact cores in BCDs 
are highly unlikely to be triggered in interacting dwarfs, if 
the underlying old stars have spherical distributions rather 
than disks. The interacting dwarf with a compact core in 
the model M2 shows an exponential radial density profile in 
its outer diffuse stellar component dominated by old stars. 
The time scale for which SFR is rather high (> O.IM© yr~^) 
in merging and interacting dwarfs is quite short (an order 
of ^ 0.1 Gyr) in the present models. Therefore, the time 
scale for merging and interacting dwarfs to be observation- 
ally identified as BCDs with very strong emission lines can 
be also short (an order of ~ 0.1 Gyr). 

The present models predict that gaseous pressure dur- 
ing starbursts triggered by galaxy interaction and merging 
can become very high (Pgas > 10^ fee) so that SCs can be 
formed in BCDs. For example, the mass fraction of SCs is 
0.22 in the model Ml, though not all of these may evolve into 
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Figure 5. The same as Fig. 2 but for the gas-poor model M4. 
Note that the central part of the merger remnant is not dominated 
by new stars in this model, which means that a BCD is highly 
unlikely to be formed in gas-poor merging even if the merger 
progenitor dwarfs have extended gas disks. 

GCs finally owing to internal and external destruction pro- 
cesses after their formation. This result clearly suggests that 
BCDs formed from dwaxf-dwarf merging can contain young 
GCs owing to tidal compression of gas during the merging. 
Furthermore, it is highly likely that nuclear star clusters in 
BCDs can be formed from nuclear starbursts during dwarf- 
dwarf merging. We thus suggest that the observed presence 
of young SCs and blue stellar nuclei in BCDs is possible 
evidence for BCD formation via dwarf-dwarf merging. 

Lastly, we describe the results of the gas-poor model M4 
in which a BCD can not be formed. Although the SFR can 
become high (~ O.6M0 yr~^) in the final phase of dwarf- 
dwarf merging, the total mass of new stars formed in the 
triggered starbursts is significantly smaller (--^ 4 x lO'^M©) 
than that of old ones (~ 6 x 10** M0). As a result of this, the 
central part of the merger remnant can not be dominated 
by new stars. As shown in Fig. 5, surface density of old stars 
within the central 1 kpc of the remnant are systematically 
much higher than that of new ones: more than an order of 
magnitude differences between old and new stars at 7? ---^ 1 
kpc. The low-mass density of new stars, combined with a 
small amount of gas mass (~ 2 x lO^M©), can not allow us 
to claim that the remnant of the gas-poor merger in this 
model is classified as a BCD. These results imply that only 
gas-rich dwarf-dwarf mergers can finally evolve into BCDs. 



4 DISCUSSION AND CONCLUSIONS 

The present study has shown that BCDs can be mergers be- 
tween dwarfs with larger fractions of gas and extended gas 
disks. The present study has also shown that (i) BCDs ap- 
pear to be compact, because the compact cores dominated 
by young stellar populations can be formed from starbursts 
triggered by merging, and (2) they also appear to be very 
blue, because the dominant stellar populations are young 
and formed from metal-poor gas during merging. Further- 
more our models have shown that extended HI gas disks of 
merger progenitor dwarfs are responsible for the formation 
of BCDs with mcissive HI envelopes surrounding blue com- 
pact cores. 

Recent HI observations have found that (i) some faint 
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galaxies have high gas mass fractions (3 < Mm /Lb < 27) 
and (ii) fainter galaxies are more likely to have higher 
Mm /Lb (e.g., Warren et al. 2006). Furthermore, the Faint 
Irregular Galaxies GMRT Survey (FIGGS) has recently re- 
vealed that the mean value of Hl-disk-to-stellar-disk ratios 
(sg) is about 2.4 for dirrs with a median Mb ~ — 13 mag 
(Begum et al. 2008). These observations suggest that galaxy 
merging that forms BCDs (i.e., mergers with high Mm/ Lb 
or /g and large Sg) can not be so rare among all types of 
dwarf-dwarf merging. 

Although some BCDs with very metal-poor, and young 
stellar populations were observationally suggested to be true 
young galaxies that are currently forming (e.g., Thuan et 
al. 1999), they have now being observed to have older stel- 
lar populations (e.g., Aloisi et al. 2007). The present study 
has shown that young stellar populations in BCDs can be 
very metal-poor, because they originate from extended gas 
disks with possibly pristine gas in merger precursor dwarfs: 
the presence of very metal-poor stellar populations in BCDs 
does not necessarily imply the very early epochs of their 
formation. 

Gil de Paz et al. (2003) investigated morphological 
types for 111 BCDs with Ha emission and revealed that 
only 10 (~ 9%) can be classified as "iI,M" (i.e., mergers). 
Our numerical results suggest that the "il" BCDs with ir- 
regular outer halos and off-center nuclei, which consists of 
35% of the BCDs (Gil de Paz et al. 2003), can be formed 
from dwarf-dwarf merging. It would be also possible that 
even "iE" BCDs with outer diffuse elliptical halos yet with- 
out clear signs of merging are merger remnants with their 
outer halos already dynamically relaxed. Although "il C" 
BCDs with cometary morphologies can be ongoing mergers 
with long tidal tails viewed from edge-on, such morpholo- 
gies can be formed from other physical processes such as 
ram pressure stripping. 

The present simulations have shown that evolution from 
BCDs formed from dwarf-dwarf merging into gas-poor dEs 
is highly unlikely owing to the presence of extended gas 
disks in the BCDs. Tajiri & Kamaya (2002) showed that 
morphological evolution from BCDs with HI envelope into 
gas-poor dEs is unlikely, because stellar feedback effects in 
BCDs are not strong enough to blow off their HI envelope. 
This observation combined with the present numerical re- 
sults therefore implies that BCDs are likely to evolve into 
gas-rich, nucleated dIrrs after significant fading of their blue 
compact cores. Papaderos et al. (1996) found significant dif- 
ferences in optical structures between BCDs and dIrrs and 
thus suggested that BCDs can evolve into dIrrs only if BCDs 
can change optical structures of their underlying stellar pop- 
ulations. By using fully consistent stellar population models 
combined with chemodynamical simulations of BCD forma- 
tion in our future papers, we discuss in detail whether the 
optical structures of BCDs formed from merging can change 
owing to fading of their young stellar populations so that the 
final structures are more similar to those of dIrrs (Papaderos 
et al. 1996). 
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